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Abstract Fix two points x,x G S 2 and two directions (without orientation) 77, 77 of the velocities in 
these points. In this paper we are interested to the problem of minimizing the cost 

Jb\ = J (»y(t)(7(t), 7(*)) + ^(t)»r(t)(7(t), 7W)) * 

along all smooth curves starting from x with direction 77 and ending in x with direction fj. Here g 
is the standard Riemannian metric on £ 2 and K 1 is the corresponding geodesic curvature. 
The interest of this problem comes from mechanics and geometry of vision. It can be formulated 
as a sub- Riemannian problem on the lens space L(4, 1). 

We compute the global solution for this problem: an interesting feature is that some optimal 
geodesies present cusps. The cut locus is a stratification with non trivial topology. 
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1 Introduction 

Fix two points x and x on a 2-D Riemannian manifold (M, g) and two directions rj G PT X M and f\ G PT X M . 
Here by PT X M we mean the projective tangent space at the point x, i.e. the tangent space T X M\ {0} with 
the identification v\ ~ ^2 if there exists a G M. \ {0} such that v\ = 0^2 . Given a vector v G T X M\ {0} we call 
Dir (i?) its direction in PT X M. We are interested in finding the path 7 : [0, T] — > M minimizing a compromize 
among length and geodesic curvature of a curve on M with fixed initial and final points x,x, and fixed initial 
and final directions of the velocity 77, f). More precisely we consider the minimization problem: 



along all smooth curves satisfying the boundary conditions 7(0) = x, Dir (7(0)) = 77, j(T) = x, Dir (j(T)) = 77, 
as in Figure [TJ Here i^ 7 (t) is the geodesic curvature of the curve 7(t), /3 is a fixed constant and the final time 
T is fixed. 



Notice that requiring 7 smooth (indeed, 7 G C 2 is enough) is necessary for the geodesic curvature being 
defined in all points. However, we will formulate the problem in a wider class to be able to apply standard 
existence results and the Pontryagin Maximum Principle (PMP in the following, see for instance [1J[2T]), that 
is a first order necessary condition for optimality. 

The problem stated above is extremely difficult in general. Indeed, one can see (see Section [2]) that on 
the projective tangent bundle PTM the minimization problem (pQ) gives rise to a contact 3-D sub-Riemannian 
problem for which global solutions are known only for few examples. The typical difficulties one meets in such 
problems are the following: 

1. One should apply first order necessary conditions for optimality (PMP) and solve an Hamiltonian system 
that generically is not integrable, to find candidate optimal trajectories (called geodesies). 

2. Even if all solutions to the PMP are found, one has to evaluate their optimality. 

In this paper we present the solution of this problem on the sphere S 2 , that we call the projective Reeds- 
Shepp car with quadratic cost. In this case the Hamiltonian system given by the PMP is Liouville integrable and 
one is faced to the problem of evaluating optimality of the geodesies. This second problem is solved with the 
computation of the cut locus K x , that is the set of points reached by more than one optimal geodesic starting 
from x. For our specific problem, indeed the cut locus coincides with the set of points in which geodesies lose 
optimality. 

The interest of this problems comes from mechanics, and more recently from problems of geometry of vision 
[HI [20l [22] . Indeed, in the spirit of the model of visual architecture due to Petitot, Citti-Sarti and Agrachev, 
the minimization problem stated above in the case M = R 2 is the optimal control problem solved by the visual 
cortex VI to reconstruct a contour that is partially hidden or corrupted in a planar image. The case M = S 2 
can be seen as a modified version of this model, taking into account the curvature of the retina an/or the 
movements of the eyes. 




(i) 




Figure 1: A curve 7 satisfying the boundary conditions. 
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1.1 Examples of problems of type (pQ) 



1.1.1 The planar case 

We present the case M = R 2 endowed with the standard Riemannian metric: it is being studied in [19]. In 
this case the problem (pQ) is equivalent to the following optimal control problem 

Jq \u\ + (3 2 u 2 )dt — > min, 




(z(0), 2/(0), 0(0)) = (*o,2/o,0o), (x(T),y(T),0(T)) = (x 1 ,y 1: 6 1 ). 

The dynamics coincide with the Reeds-Shepp car (see [22]). except for the fact that here G R/7T. For this 
reason we call the problem (pQ) the projective Reeds-Shepp car on a Riemannian manifold, with quadratic 
cost. 



1.1.2 The case M = S 2 

The case M = S 2 , endowed with the standard Riemannian metric, can be seen as the compactified version of 
the problem on the plane. It is the optimal control problem 

Io( u i + P 2 U2)dt — > min, 



(a(0),&(0),£(0)) = (a ,6o,&), (a(T), 6(r),$(T)) = (ai,&i,&), 

where (a, 6) are standard spherical coordinates on the sphere (see ((!])) and £ = Dir (7) . Also in this case the 
Hamiltonian system given by the Pontryagin Maximum Principle is Liouville integrable and one is faced to the 
problem of evaluating optimality of the geodesies. 

This problem has special features. Indeed, PTS 2 is the 3-D manifold called lens space L(4, 1), that can be 
seen as a suitable quotient of SU (2) by a discrete group. Moreover, when (3 = 1 this problem is the projection 
of a left- invariant sub- Riemannian problem on SU(2), called k © p problem. For SU(2) an explicit expression 
for geodesies is given and we computed the cut locus and the Carnot-Caratheodory distance in [9] (results are 
recalled in Section IcTT]) . As a consequence, one can easily find geodesies for the problem on PTS 2 as projections 
of the geodesies for SU(2). An interesting feature is that the projections of these geodesies on S 2 present cusps 
in some cases, see Figure [2] Observe that in cusp points the tangent direction is well defined. 







(•; 


= ui 







cos(e) \ 

sin(g) 
sin(a) 
■ cot(a) sin(£) J 





Figure 2: A geodesic with a cusp. 
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L(4, 1) can be described tope-logically as follows: consider a full 3-D ball B := 

and define the equivalence relation w on it as follows: ~ is reflexive; the points (x^x^x^ 
dB n {a; 3 > 0} and (x^ , x 2 , x 3 ) G 95 — dB D {^3 < 0} are identified when 

4 = -3-and ( X x i) = (° 1 ~~o){ X x;)> 
as in Figure 03 The manifold B/~ is topologically equivalent to L(4, 1). 




Figure 3: Identification rule on dB. 



Observe that there exist different identification rules on SU (2) such that the quotient manifold is topolog- 
ically L(4, 1) and the induced sub-Riemannian structure is well defined: indeed, in [9l Section 4] we already 
endowed L(4, 1) with a sub-Riemannian structure that is different from the one we present in this paper. The 
structure chosen in this paper is the unique for which the minimal length problem on L(4, 1) coincide with 
the problem (pQ). As a consequence, also the cut loci are different in the two cases. See Remark [3] for further 
details. 

The main result of this paper is the computation of the cut locus if [id], presented in Theorem [H It is a 
stratification and it has nontrivial topology. Indeed, it is the union of two sets K [Id] = K loc U K a ^ m , described 
as follows: K loc is topologically equivalent to S 1 without the starting point [Id]. K sym is given by the union 
of two 2-D manifolds glued along their boundaries 

K sym =K a uK b 



where 

• The manifold K a is contained in the set |[#] | g = ^ °^ ^_ ^ with Re (/?) = o| and it is home- 

omorphic to a 2-D disc with only two points identified on the boundary. Observe that the boundary is 
thus homeomorphic to two circles S* 1 glued in a point (we call them and s 1 ). A part of K loc (namely, 
the farthest part from the starting point) is contained in K a . 

• The manifold K b is topologically equivalent to a 2-D disc whose boundary S* 1 is glued to the boundary 
of K a as presented in figure [H 

A representation of the topology of the cut locus is given in figure [U 
We give a picture of the cut locus in Figure El Here: 
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Figure 4: A scheme of the topology of the cut locus: points Pi are identified, boundary of disks are identified 
according to arrows. 



• The 1-D stratum K loc is represented by two semi-diameters (one is the left-right line on the northern 
hemisphere, the other is the front-rear line on the southern hemisphere) without the poles (representing 
[Id]): due to identification rule ~, the two lines are identified. 

• The 2-D stratum K a is represented by four "triangles" on the boundary of the sphere. The two dark gray 
triangles (on the front) are identified . Similarly, the two light gray triangles (on the rear) are identified. 

• The 2-D stratum K b is not represented for a better comprehension: it can be seen as a surface inside the 
sphere whose boundary is the closed line given by the concatenation of segments SlZSIZsrsr. 

The paper is organized as follows: in Section [2] we define the mechanical problem on an orientable Rieman- 
nian manifold M and write the optimal control problem in the case M = S 2 . In Section [3] we define L(4, 1) and 
a sub-Riemannian structure on it such that the minimal length problem on L(4, 1) coincide with the mechanical 
problem stated above. Finally, in Section [4] we give the solution of the problem. 

2 The problem of minimizing length and curvature on a 2-D Rie- 
mannian manifold 

In this section when we speak of directions we always mean directions without orientations. 
Let M be a 2-D manifold smooth manifold and PTM its bundle of directions, i.e., 

PTM := \J {PT X M}, where 

xeM 

PT X M := {(T X M \ {0})/ - where v± ~ v 2 if and only if 3 a e R \ {0} | v 1 = av 2 } • 

Assume now that M is orientable and Riemannian, with Riemannian metric g and let (xi,^) be a local 
orthogonal chart on an open set U C M, i.e. a chart such that, 

E( Xl ,x 2 )(d Xl ,d X2 ) = on U. 

In the following the symbol || • || indicates the norm with respect to g( Xl , X2 )- 
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In this case there is a canonical way to construct a chart (x±, £2, £)> w ^ n £ ^ W 71 " on PTU , as follows: write 



17 = 



■ ^21 



l& a |, 

and define £ as the angle between i? and d Xl /\\d Xl ||, i.e. £ := arctan (j^j • 

In the following we call the chart (xi,£2,0 a lifted- orthogonal chart. Notice that \\d Xl \\ = 1 and H^ll = 1 
on U, implies g flat on U. 

Given a smooth curve 7(.) : [0,T] -> M, its Z£/fc in PTM is the curve (7(.), £(■))> where G PT 7(£) M is 
the direction of 7(t). Writing 

<y(t) = 71 (t)0 Xl +72(i)^ 2 = (71 WH^ 11)^/11^ || + (Ut)\\d X2 \\)d X2 /\\d X2 \\, 



then 

Given a smooth curve T(.) : [0,T] — > PTM, we say that it is horizontal if it is the lift of a curve on M. If 
r(-) = (7i(-)>72(-)> £(■)) are its components in a lifted orthogonal chart, then T(.) is horizontal if and only if 
condition ([2]) holds. 

The requirement that a smooth curve T(.) in PTM is horizontal, is equivalent to require that its velocity 
belongs to a 2-D distribution (see Appendix A). Let us describe this fact in detail: in a lifted orthogonal chart 
we have T(t) = ji(t)d Xl + 72(^)^2 + ano ^ condition ([2j) is equivalent to 

/ cos(C) \ 



f(t) eA(T(t)), where A(xi,x 2 ,0 := span{*i, F 2 } , Fj = ' 



11% II 
sm(Q 

113*, II 



F, = 



\ Z(x 1 ,x 2 ,0 j 
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The function E(xi,x 2 ,^) is an arbitrary smooth function. We choose H such that the vector field F\ on PTM 
is the geodesic spray, i.e. the infinitesimal generator of the geodesic flow (see [HI [25]). 

Notice that span{Fi(xi, x 2 , £), F 2 (x x , x 2 , ,[Fi,F 2 ](xi, x 2 , £)} = T (x u x 2 £)( PTM )i nence (PTM, A) is a 
3-D contact distribution. 



2.1 A sub-Riemannian problem on PTM 

There is a natural sub-Riemannian problem associated with (PTM, A). It comes from the problem of mini- 
mizing a compromise among length and geodesic curvature of a curve on M with fixed initial and final points, 
and fixed initial and final directions of the velocities. More precisely we consider the following problem: 

(P) Fix x, x e M, r] e PT X M and f) e PT X M . Minimize 

m = £ v /g 7( t)(7W,7W)+^ 2 ^ (£) g 7 (t)(7W,7W) dt (3) 

along all horizontal smooth curves T(.) = (7O) >£(■)) : [0,T] —> PTM satisfying the boundary conditions 
T(0) = (x,rj), T(T) = (x,f)), i.e. 7(0) = x, Dir(7(0)) = 77, 7 (T) = x, Dir (<y(T)) = fj. Here K l{t) is the 
geodesic curvature of the curve j(t), and (3 is a fixed nonvanishing constant. 

Problem (P) is illustrated in Figure [TJ Its equivalence with the problem (pQ) (in which the square root 
is absent) is eplained below. The cost ([3]) is the most natural cost associated with (PTM, A). Indeed it is 
invariant by change of coordinates and by reparameterization of the curve 7. Moreover, as it will be clear from 
the example below, it can be interpreted as a length in PTM. 

The term g 7 (t) (7(^)5 j(t)) corresponds to the Riemannian length on M, while the term /^ 2 if 2 ^g 7 ( £ )(7(t), j(t)) 
corresponds to the geodesic curvature. The presence of g 7 (t)(7(t), j(t)) in this second term guarantees the in- 
variance by reparametrization. The constant f3 fixes the relative weight. It may happen that there is a natural 
choice for (3 or even that (3 is not relevant (see an example on R 2 below), depending on the specific problem. 

One can check that g 7 (t)(7(t), j(t)) + P 2 K 2 ^g^^(j(t), j(t)) defines a positive definite quadratic form on 
A, i.e. it endows (PTM, A) with a sub-Riemannian structure (see Appendix [Aj) . Indeed, if x = (x\,x 2 ,£) is 
a lifted orthogonal chart and H is such that F\ is the geodesic spray, we have that problem (P) becomes the 
optimal control problem: 



/ cos(0 \ 

U1F1 + u 2 F 2 , Ft = 



11% II 
sin(C) 



\\d X2 \\ 
\ S(xi,x 2 ,0 J 




m = f 

Jo 



P 2 Uo dt — > min 



that is the minimal length problem in the sub-Riemannian manifold PTM. 

Consider now the problem of minimization of the energy functional E[T] := J Q T u\ + f3 2 u\ dt with fixed 
final time T and fixed starting and ending points. If T is a minimizer of E, then it is a minimizer of J and 
its lifted velocity v(t) := \j u\ + (3 2 u\ is constant. On the other side, a minimizer r of J parametrized with 
constant lifted velocity v is also a minimizer for E with T = J\T]/v. For details, see Appendix lAl 

Finding a complete solution for the problem (P), i.e. finding the optimal trajectory connecting each initial 
and final condition, is a problem of optimal synthesis in dimension 3 that is extremely difficult in general. The 
case M = R 2 with the standard euclidean structure is being studied in [19]. The complete solution in the 
case M = S 2 with the standard Riemannian metric is given as the solution of the sub-Riemannian problem on 
L(4, 1) presented in Section 03 



2.1.1 The planar case 

Consider M = R 2 with the standard euclidean structure. In this case PTR 2 ~ R 2 x S 1 and a lifted orthogonal 
chart is (x\,x 2 ,£), where (x\,x 2 ) are euclidean coordinates on R 2 . In this case, since \\d Xl \\ = \\d X2 \\ = 1, 
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writing j(t) = (xi(t),X2(t)), we have 



5(7? 7) = #1 + %2-> an d ^7 



-X2 X\ + Xi X2 



, where = arctan — 



Xl (t)J ' 



Hence J\F] = L y x 2 + x?, + ft 2 !; 2 dt. We choose #2, = 0, then Fi is the geodesic spray and we have 



/ COS(0 \ / \ f T 

A( Xl ,x 2 ,0 :=span{Fi,F 2 }, where F 1 = I sin(£) I, F 2 = I I , JfT] = ^ y'u? . 



[3 2 ul dt. 



The associated optimal control problem is 



x\ = ui cos(£) T 

x 2 = u\ sin(£) m 1 ,« 2 €R, / y u i+ P 2u \ dt — > min . 

This problem can be seen as a left-invariant sub-Riemannian problem on the group SE(2)/ ~, the group of 
rototranslations of the plane 

cos(a) — sin(a) x\ 

SE(2) := ^ [ sin(a) cos(a) X2 
1 



a G JR/27T, aji G 



endowed with the identification rule (a, xi, #2) ~ (a + 7r, #1, #2)- In this case the constant /? can be set to 1, 
since it becomes irrelevant after the transformation (xi,X2) — > (/3xi,(3x2). This problem is being studied in 

eg. 

2.2 The case M = S 2 

Let S 2 = {(2/1,2/252/3) G M 3 I 2/1+2/2+2/3 = 1} an( ^ ^ ( a 5 fr) ^ e ^ ne spherical coordinates (see Figure [6|): 



2/i = sin(a) cos(6) 

2/2 = sin(a) sin (6) a G [0, 7r], 6 G [0, 2+). 
2/3 = cos(a), 



(4) 





y 3 










^/^+ii^ 





Figure 6: Spherical coordinates on S 2 . 
In these coordinates the standard Riemannian metric on S 2 (induced by its embedding in R 3 ) has the form 



sJ(a) 



The geodesic curvature of a curve j(t) = (a(£), b(t)) is: 



y sin(a) 2 ^2 cot(a) d 2 b + cos(a) sin(a) b 3 — ba + ab 



cos(a) b 



^d 2 + sin(a) 2 6 2 ^ 



d 2 + sin(a) 2 6 2 A/d 2 + sin(a)^ 6 2 



\2 



with £ = arctan 



/ 6 sin(a) \ 

V * / ' 



Hence J[r] = J Q J d 2 + sin (a) 6 2 + (3 2 (cos(a)b + £) 2 . In this case Fi is the geodesic spray choosing 
S(a, 6, £) = — cot(a) sin(£) for which we have 



A(a, &,£):= span {Fi, F2} , where Fi 



/ cos(0 \ 
sin(Q 
sin(a) 
y — cot(a) sin(£) J 



f 2 = ^ j , j[r] = jf ^7^5 a. 



The associated optimal control problem is 



' a = u\ cos(£) 

• sin(f) 
= u\- 



sin(a) 

£ = — u\ cot(a) sin(£) + 1x2, 



2 dt — > min . 



(5) 



Remark 1. Observe that this optimal control problem admits trajectories for which d = 6 = 0. These trajec- 
tories in PTS 2 represent the possibility of changing the direction of the tangent vector on the sphere S 2 on a 
fixed point; mechanically, it represents the possibility of rotation on itself. 



2.2.1 Existence of minimizers 

To discuss the problem of existence of minimizers for the problem on the sphere S 2 , let us come back to the 
original problem of minimizing the cost (pQ), where M = S 2 with the standard Riemannian metric, along all 
smooth curves 7 : [0,T] — > S 2 . Here the class of smooth curves has been chosen to give a meaning to the 
geodesic curvature in the whole interval [0,T] (indeed C 2 ([0,T]) was enough). 

Unfortunately the class of smooth (or C 2 ) curves is too small to apply standard existence theorems for 
minimizers and first order necessary conditions for optimality (the PMP). We have to deal with a larger class 
of function. 

This class arises naturally after the lift in PTS 2 , where the minimization problem takes the form ([5]). Here 
it is natural to look for minimizers (a (.),&(.),£(.)) that are absolutely continuous and corresponding to controls 
u\(.) , U2(.) belonging to F 2 ([0,T]). In this class, standard existence theorems can be applied to guarantee 
existence of minimizers (see for instance [24j Theorem 5.1]). 

However the PMP works in the smaller class of L°° controls (corresponding to Lipschitz curves (a(.), &(.), £(.))); 
hence, before applying PMP, one needs to prove that minimizers belong to this smaller class. 

Thanks to a theorem of Sarychev and Torres (see [24j Theorem 5.2]), one can prove that for the problem ([5j) 
all minimizers satisfy the conditions given by the PMP and that normal ones correspond to controls belonging 
to L°°([0, T]). Since in our case there are no abnormal extremals (the problem (j5j) is 3-D contact), it follows 
that all optimal controls are indeed L°°([0,T]). 

Finally, all trajectories satisfying the PMP are solutions of an Hamiltonian system, that in our case is 
analytic. Hence all minimizers are analytic and therefore smooth. 

We conclude that it is equivalent to solve the original problem in the class of smooth curves or the lifted 
problem in the class of L 2 controls. A similar treatment has been presented with all details in [24] for a similar 
problem but in the presence of a drift. 
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3 A sub-Riemannian problem on L(4, 1) 



In this section we define a k p sub-Riemannian structure on SU (2) and compute its cut locus and the sub- 
Riemannian distance. For details on Sub-Riemannian geometry and k®p manifolds, see Appendix [Al We then 
define a sub-Riemannian structure on L(4, 1) induced by the one on SU (2) and prove that the sub-Riemannian 
problem on L(4, 1) coincide with the optimal control problem ([5|) on PTS 2 . 

3.1 The k © p problem on SU(2) 

The Lie group SU (2) is the group of unitary unimodular 2x2 complex matrices 

SU(2) = {( _^ g ) G Mat(2,C) | \a\ 2 + \f3\ 2 = 

It is compact and simply connected. The Lie algebra of SU (2) is the algebra of antihermitian traceless 2x2 
complex matrices 

su (2) = | ( 1^ -ia) e Mat ( 2 > C ) ' a E M ' & E C } ' 
A basis of su(2) is {pi,P2, k} where 



1/0 1 \ 1 / i \ , 1 / i 

P1 = H-1 J P2= 2{i J fc= H0 -< 

whose commutation relations are [pi,P2] = k [p2,k] = pi [k,pi] = P2- Recall that for su(2) we have 
Kil(X, Y) = 4Tt(XY) and, in particular, Kil(^,pj) = —2Sij. The choice of the subspaces 

k = span {k} p = span {pi , p2 } 

provides a Cartan decomposition for su(2). Moreover, {^1,^2} is an orthonormal frame for the inner product 
< •, • >= — ^Kil(-, •) restricted to p. 

Defining A (g) = gp and g g (vi, V2) =< g~ x v\,g~ x V2 >, we have that (SU(2), A, g) is a k®p sub-Riemannian 
manifold (for details see Appendix IA.3|) . The sub-Riemannian manifold SU(2) is thus endowed with the 
standard definition of sub-Riemannian length and distance. 

3.1.1 Expression of geodesies 

k p manifolds have very special properties: there are no strict abnormal minimizers, hence all the geodesies 
starting from a point go are parametrized by the initial covector. The explicit expression of a geodesic starting 
from go is 

g(t) = goe( Ak + A ^e- Akt . (6) 

with Ak G k, A p G p. The geodesic is parameterized by arclength when < A pi A p >= 1. This condition defines 
a cylinder C T* Q M. 

For the SU(2) case, we compute the geodesies starting from the identity using formula (|6]). Consider an 
initial covector A = A((9, c) = cos(0)pi + sin(f9)p2 + ck G A^. The corresponding exponential map is 

Exp(<9, c, t) := Exp(A(0, c), t) = e ^os(e) Pl+ sin(e) P2+ ck)t e -ckt = ( 

with 

a 



(3 a 

(f )sin(Vl+c 2 |) ^(ct\ /1 1 J2t\ 1 • f ccos(f )sin(Vl+c 2 |) 



Vl+c 



+ cos(f ) cos(VTT^I) + i ^-(f^f+^ il _ sin (|) cos(VTT^I)) 



P J \ Mf + g)+<sin(f + 0)) 

Another property of k © p manifolds is the existence of a solution of the minimal length problem (see 
Remark [T2]) , i.e. for each pair g,h G M there exists a trajectory 7 steering g to h and minimizing the length: 
= d(g, h). Recall that this minimizing trajectory is a geodesic. 
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3.1.2 The cut locus and distance for SU{2) 

In this section we recall the formula of the cut locus for the k © p manifold SU (2) and of the sub-Riemannian 
distance. Both the results are proved in [9]. 

Theorem 1. The cut locus for the k ® p problem on 577(2) is 

K ld = e k \Id = {e ck | c G (0,4tt)}. 
The cut locus is topologically equivalent to a circle S 1 without a point, the starting point Id. 

Theorem 2. Let g = ( °^ - ) ^ SU(2). Consider the sub-Riemannian distance from Id defined by the 



-p a 

k p structure on SU (2). It holds 



dQd,g) = l VargW(2.-arg(.)) z/ = 

where arg (a) G [0, 27r] and ^(a) = ^ where t is the unique solution of the system 

-f + arctan (^j^y tan (^^)) = arg(a) 



Remark 2. Notice that V a, ft, /fe G C,|ft| = |ft | we have d [id, y ^ J J = d [id, y ^ J J = d ^Id, ^ ^ 
3.2 A sub-Riemannian problem on L(4, 1) 

We define the 3-D sub-Riemannian manifold L(4, 1) as a quotient of the k © p manifold SU (2) defined above. 
L(4, 1) is a lens space: for more details about lens spaces L(p, q) see [23]. We prove that the quotient and the 
sub-Riemannian structure are compatible. We then compute the distance d(., .) on L(4, 1). 

3.2.1 Definition of L(4,l) 

We define coordinates (a, 6, c) on SU(2) as follows: 
Proposition 1. £ac/i g G SU(2) can be written as 

g = e -b P 2 e ap le cp 2 ^ 

with a G [0, 7r]j b G [0, 27r), c G R/(47r). The value of a is uniquely determined by g. 
If a ^ 0, a 7T, then the values of b and c are uniquely determined by g. 
If a = ; then the value of c — b mod 4tt is uniquely determined by g. 
If a = 7T ; £/ien £fte mfo/e of c-\-b mod 47r zs uniquely determined by g. 

Proof. Observe that e -^ 2 e opi e cp2 = ^ _^ ^ ^ with 

fa = cos (f ) cos + z sin (f ) sin (^) 

\/3 = sin (f ) cos (^) + i cos (§) sin (^) . 

The result follows from direct computation. □ 
In the following we will write g = (a,b,c) for g = e~ bp2 e api e cp2 , even though the vector (a, 6, c) is not 
uniquely determined by g. 
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Proposition 2. Define the following equivalence relation on SU{2): g\ w g^ if there exist (ai, &i, ci), (<22, 62, C2) 

{ai = a 2 
61 = 62 
Ci = C2 mod 7T. 

T/ie 3-1} manifold SU(2)/ w zs £fte /ens .space L(4, 1). 

Before proving the proposition, we recall the standard definition of L(p, g) with p, g G Z\ {0} coprime. Let 

S 3 = {(x 1; . 2 )eC 2 I H 2 + |x 2 | 2 = i} 

and define on S* 3 the following equivalence relation: (#1,2:2) ~ (2/1? 2/2) if there exists wGC p-th root of unity 
(i.e. cj p = 1) such that 

x 2 y V w V ^ y 2 J ' 
The quotient manifold S 3 / ~ is the lens space L(p, g). 

Sf7(2) -> S 3 , 
Proof of Proposition^ Consider the isomorphism : f a (3 



-(3 a 



( a+/3 Re(cQ+iIm(/3)-Re(/3)-iIm(oQ \ 
I V2 ' 72 ; 



A straightforward computation gives that #i, #2 £ £^ (2) are equivalent (#i « #2) if and only if 0(#i), ^(^2) £ £ 3 
are equivalent ((j)(gi) ~ 4>{d2)) with respect to the equivalence relation defined by p = 4, g = 1. Hence the 
manifolds SU(2)/~ and 5' 3 /^ are isomorphic. Thus SU(2)/~ is a 3-D manifold, the lens space L(4, 1). □ 

3.2.2 Sub-Riemannian structure on L(4, 1) 

We now define a sub-Riemannian structure on L(4, 1) induced by the quotient map. 

Proposition 3. The sub-Riemannian structure on SU(2) given in \3.1\ induces a 2-dim sub-Riemannian struc- 
ture on L(4, 1) via the quotient map 

n . SU(2) - L(4,l), 

9 ^ \g], 

i.e. 

• the map 

A : [g] i-> II* (A(ft)) C T b] L(4, 1) wz*A A e [#] 
25 a 2-dim smooth distribution on L(4, 1) £fta£ zs Lie bracket generating; 

• g[s](^M^*) = : = ( v i w )h w tth h e [#]> G T h SU(2) : II* (v) = v*, II* (w) = w* is a smooth 
positive definite scalar product on A. 

Proof. The role of the map II and 11*^ is illustrated in the following diagram 

T g SU(2)^^T [g] L(4,l) 



577(2) 5-^L(4,l). 

The map II is a local diffeomorphism, thus II* |g : T g SU(2) — > T^L(4, 1) is a linear isomorphism, hence 
II*! (A(g)) is a 2-dim subspace of T^L(4, 1). 
The following statements: 

• the distribution A([g]) is well defined, i.e. V hi, /12 £ [#] we have 

n, |fc (A(h 1 ))=n* lh (A(h 2 )); 
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• the positive definite scalar product (v*, w*)[ g ] is well denned, i.e. V fti,ft2G[g], vi,wi G Tf ll SU(2), V2,W2 G 
T h2 SU(2) such that 11*^ (vi) = 11* j (v 2 ) and 11*^ (wi) = 11*^ (w 2 ) we have (vi,wi) hi = (v 2 ,w 2 ) h2 ; 

are consequences of the following lemma: 

Lemma 1. Let fti, ft 2 £ [#] wztft fti = (a, c) and ft 2 = (a, c + kir), k G Z. Tfte map 

* P ~> P- 

ra = raipi + ra 2 p2 >-> (-l) fe raipi +ra 2 p2 

zs bijective and it is an isometry w.r.t. the positive definite scalar product ( , ). Tfte following equality holds 
V m G p 

Proof. A direct computation gives that (j) is bijective and that it is an isometry. 
Consider now h\e trn and ft 2 e* n with m = m\p\ + m2P2 G p, n = n\p\ + 712^2 G p. We have 

hi e tm = e -b P2 e a P2 e c P2 e tm = ( a(a, 6, c, t, m) /3(a, 6, c, t, m) 

-P(a,b,c,t,m) a(a,b,c,t,m) 



with 



a(a, 6, c, t, m) = -Ko(£j m)a(a, 6, c, 0, 0) + (—mi + irri2)Ki(t, m)/3(a, 6, c, 0, 0), 
(3(a,b,c,t,m) = i^o(^ 5 m)/3(a, b, c, 0, 0) + (mi + ira 2 )-Ki(t, m)/3(a, fr, c, 0, 0), 



K (t,m) := cos £ , K^m) := 



sin 



2 / \/ m l + m 2 

Observe that ft 2 e tn = e ~ bp2 e ap2 e cp2 + k7TP2 e tn ~ f ^ 2 ^ 2 ) with 

«2 = K (t,n)a(a,b,c,0,0) + ((-l) fe+1 ni + in 2 )Ki(t, n)p(a, 6, c, 0, 0), 
/3 2 = #o(*, rc)/?(a, 6, c, 0, 0) + ((-l) fc ni + m 2 )Xi(t, n)/?(a, 6, c, 0, 0), 



hence [ftie^ m ] = [ft2e £n ] if both the following conditions are satisfied: 



' K (t,m) = K (t,n) 

-mi + im2 = (— l) fc+1 ni + m 2 
jri\ + im 2 (— l) fc ni + m 2 . 
These are verified when n = 0(m). Thus ^ [ftie £m ] = [ft 2 e*^ (m) ] • □ 

Since II is a local diffeomorphism, V g G SU(2) 3 B(g) such that the map ^ : T B ^SU{2) — > 

Tg([p])L(4, 1) is a diffeomorphism, thus A is smooth, Lie bracket generating, and (v*,^*)^] is smooth as a 
function of [g]. □ 

Remark 3. Observe that other definitions of a sub-Riemannian structure on L(4, 1) induced by the one on 
SU(2) are possible: for example, in [9] we defined the following identification rule on SU{2): 

% ^ ) ~ ( ) if exists co e C such that uo p = 1 and 1^/1 = 1^ ^ I f °> I (8) 

-(5 a J \ -f5' a' J \ f3' J \ u q J \ (3 J 

and observed that the sub-Riemannian structure defined on SU(2) in 13.11 induces a sub-Riemannian structure 
on the manifold L(p,q) := SU(2)/ <~. 

The sub-Riemannian structures are locally isomorphic, but globally the structure may vary: indeed, in the 
following we will prove that the cut locus for the structure defined in [3] is globally different from the one defined 
by ®. 
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Remark 4. Notice that in this case the push-forward of a left invariant vector field is not a vector field, because 
the projections (gp), 11*^ (ftp) from different points g, ft such that [g] = [ft] do not coincide. Nevertheless, 
the projections of the whole distribution II* j (A(g)), 11*^ (A(h)) coincide; then the sub-Riemannian structure 
on L(4, 1) is well defined even though the projections of the left invariant field is not well defined. 

We have standard definitions of sub-Riemannian length and distance on L(4, 1), as presented in Appendix 
[A| see ([TO]) and ([TT]) . We indicate them respectively with / and d. 

Observe that the geodesies for the sub-Riemannian manifold L(4, 1) are the projections of the geodesies for 
SU(2), due to the fact that the projection II is a local isometry. In general we have the following 

Proposition 4. Let T : [0, T] — > SU(2) be a smooth curve. Define its projection T_ : [0,T] — > L(4, 1) by 
T(t) := H(T(t)). Then T is a smooth curve and its length coincide with the length ofT: 

l(T) = l(T). 

Recall now the definition of the lift of a curve in our case and a standard result about its optimality: 

Definition 1. Let 7 : [0, T] — > L(4, 1) be a smooth curve in L(4, 1) with 7(0) = [g\. Fix g\ E SU(2) such that 
gi G [g\- The lift of 7 starting from g\ is the unique smooth curve 7 : [0,T] — > SU(2) satisfying 7(0) = g\ and 
n(7(t))= 7 (t) for anyte[0,T}. 

The length of the lift 7 in SU(2) coincide with the length of the curve 7 in L(4, 1); 

m = On)- 

Proposition 5. Let 7 : [0,T] — > L(4, 1) be an optimal trajectory steering [g] to [ft]. Fix g\ G SU(2) with 
9i £ [g]- Then the lift 7 : [0,T] — > SU(2) starting from g\ is an optimal trajectory steering g\ to j(T) G [ft]. 

We end this section giving an expression to compute the sub-Riemannian distance on L(4, 1). 

Proposition 6. The sub-Riemannian distance d on L(4, 1) can be computed via the sub-Riemannian distance 
d on SU (2) as follows: 

d([g], [ft]) = min hie[h] {d(g u hi)} 

independently on the choice of g\ G [g] . 

Proof. Consider a trajectory 7 in L(4, 1) that is solution of the minimal length problem between [g] and 
[ft]. Define its lift 7(.) starting from a fixed g\ G [g]. It steers g\ to a given h\ G [ft] and it is optimal, thus 
d([g],[h]) = d(g 1 ,h 1 ). 

We prove now that for any hi G [ft] we have d(g\, h\) < d(gi, ft2). By contradiction, assume that 3 G [ft] 
such that d(gi, I12) < d(gi, fti), thus there exists a smooth curve T : [0, T] — > SU (2) steering g\ to ft2 with l(T) < 
d(gi,hi). Consider its projection T: it steers [g] to [ft] and it satisfies l(T) < d([g], [ft]). Contradiction. □ 



3.3 Isometry between L(4, 1) and PTS' 2 

This section is fundamental in the treatment of the mechanical problem presented in Section [2j we state that 
the problem ([5]) on PTS 2 is the problem of minimal length on the sub-Riemannian manifold L(4, 1). 
We start writing the expression in coordinates of the projection S 3 — > 5 2 , the well-known Hopf map. 

Proposition 7. Consider the Hopf map 

SU(2) ~S 3 -> S 2 , 

2 (Re (a) Re (/?) + Im (a) Im (^)) \ 
2 (Re (a) Im (a) - Re \(3) Im (^)) 
Re (a) 2 + Im (/3) 2 - Re (f3) 2 - Im (a) 2 . / 

Consider the coordinates (a, 6, c) on SU(2) defined in (0) and £fte spherical coordinates (a, 6) on 5 2 defined in 
(0); £ften £fte i^op/ map expression in coordinates is 



a (3 
—/3 a 



14 



SU(2) S 2 , 
(a, 6, c) i— > (a, 6). 

Remark 5. Observe that both the (a, 6, c) coordinates on SU(2) and the spherical coordinates (a, 6) are not 
well defined in a = and a = n. Nevertheless, the map <fi is well defined and has this expression in coordinates 
also in the degenerate cases. 

The following theorem is the central result of this section: here we prove that our mechanical problem ([5j) 
coincide with the sub-Riemannian problem on L(4, 1). 

Theorem 3. The lens space L(4, 1) is diffeomorphic to PTS 2 , the bundle of directions of the sphere. The 
coordinates (a, 6, c) on £(4,1) coincide with the coordinates (a, 6, — c) on PTS 2 . The diffeomorphism maps 
horizontal curves on L(4, 1) to horizontal curves on PTS 2 and vice versa. 

The sub-Riemannian length of an horizontal curve (a(. ),&(.), c(.)) in L(4, 1) zs egim/ to t/ie cost i7[r] of the 
curve = (a(.), 6(0. - c (0) ^ PTS 2 . 

Proof. Consider any horizontal smooth curve T : (—£,£) — > SU(2) with T(0) = (a, 6, c). We have 
T(t) = e~ bp2 e api e cp2 e tv ^°^ with ?/Gp due to the fact that T is horizontal. Assume 77 = 771P1 -\-1j2P2 satisfying 
771 7^ 0. The reason of this condition will be clear in the following. 

Consider its projection 7 = ip o T : (— £, e) — > 5 2 : we have by computation 

(sin (a) cos(6) \ / cos(a) cos(6) cos(c) + sin (6) sin(c) \ 

sin(a) sin(6) +trji cos(a) sin(6) cos(c) — cos(6) sin(c) +o(t). 
cos(a) J \ — sin(a) cos(c) J 

Due to condition 771 7^ 0, 7 is smooth and differentiate in 0; the direction of its tangent vector in 
depends only on T(0) = (a, b, c). Its lift is (a(. ),&(.), £(.)) with £(t) = arccos ( ^ S1 ^ a ^ : an explicit computation 

SU(2) —> PTS 2 

gives £(t) = — arctan(tanfc)). Thus the function f : , 7 y { , 7 ' /L , UN maps horizontal 

to sw v w/ (a,6,c) 1 — ► (a,6,-arctan(tan(c))) 

smooth curves of SU(2) with 771 7^ to horizontal smooth curves of PTS 2 . 

The case 771 = can be studied as a limit case: in this case a = b = 0, thus / maps horizontal curves 
of 577(2) to the curves of PTS 2 representing rotations on a point (see Remark [T]). In this sense / maps all 
horizontal smooth curves of SU(2) to horizontal smooth curves of PTS 2 . 

Consider now L(4, 1) = SU(2)/&: the map / pass to the quotient and we have the explicit expression 

L(4,l) -> PTS 2 , 
[a, 6, c] 1 y (a, 6, — c). 

Observe that /~ is a diffeomorphism. Moreover, it maps horizontal smooth curves of L(4, 1) to horizon- 
tal smooth curves of PTS 2 . Then the differential df~ is a linear isomorphism df~ : A([g]) C L(4, 1) — > 
A(a, 6, — c) C PTS 2 that is explicitly df~(rjiU* lg (gpi) + rj2U^ (gp 2 )) = rjiFi(a,b,—c) - rj 2 F 2 (a, 6, -c) with 
g = (a, 6, c) G *Si7 (2). Then we have that the length of a curve 7 : [0, T] — > L(4, 1) satisfies 



f(a(.),6(.),c(.)) = V^+^ = '[«(■)>*(■)]■ 



□ 



4 Proof of the main result 

We give now the explicit computation of the cut locus for the problem (P) on S 2 . For this specific problem the 
first conjugate locus (see Appendix IA.2j) is completely contained in the cut locus since this is the case for the 
lifted problem on SU(2) (see [9j Sections 3.1.3, 5.1]). Hence the cut locus is the set of points where geodesies 
lose optimality (see Remark [9]). 
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Theorem 4. The cut locus if [id] f or the sub-Riemannian problem on L(4, 1) starting from the identity is a 
stratification K [ld] = K loc U if sym et/i 

if' oc = [ifid]\[Id] = [e k ]\[Id], 
K sym = G £^ | fl g * £ 3 Xi, x 2 G [Id] , a:i 7^ x 2 sizcA £/ia£ 

Vx* e [Id] xi) = d{g\x 2 ) < d{g\x*)} . 

where .) zs £/ie sub-Riemannian distance on SU(2) given in Theorem^ 

Proof. We first prove that K loc lies in the cut locus. Consider [g] G K loc : due to the definition of 
[#], there exists #1 G [#] such that #1 = e nk = ^ ^ ^ ^ with |a| = 1. Hence [g] = {#1, -#i, #2, -#2} 
/ \ 

with g 2 = I • )• Recall that d([Id], [g]) = min^*^] d(Id,g*). Applying Theorem [21 we have d(Id,gi) = 



~ / 7 \ 77: t/t n \ I 2a/ (arg(a) + 7r)(7r — arg(a)) ifaG[0,7r] 7/T , x 7/T , x 

2 Varg (a 2tt - arg (a , d Id, -</i = { /r^H Jr , <Id, <? 2 ) = <Id, -<? 2 ) = 

I 2 v / (arg(a) - 7r)(37r - arg(a)) if a e [tt, 2k) 

7r. A straightforward computation gives that the minimum is attained for d(Id, #i) or d(Id, —g\). We assume 
without loss of generality that we have d([Id], [g]) = d(Id,gi). Recall that g\ G if id, thus there exist 71,72 
two different optimal trajectories on SU(2) steering Id to g\. Thus their projections 71,72 are trajectories on 
L(4, 1) steering [Id] to [g] and satisfying l(ji) = = d(Id, #i) = d([Id], [g]) for i = 1, 2; hence they are both 
optimal. Observe that 71, 72 are distinct in a neighborhood of Id, thus 71, 72 are distinct in a neighborhood of 
[Id]. Thus [g] is a cut point, due to the fact that it is reached by two different optimal trajectories 71,72- 

We prove now that K syrn lies in the cut locus. Consider a point [g] G K sym and a fixed g* G [g\. Recall 
that d(g*,xi) = d(g*,x 2 ). Due to the existence of a solution to the minimal length problem, there exist I\ 
and T2 geodesies in SU (2) with the following property: I\ steers X{ G [Id] to g* G [g] and it is optimal. 
Observe that Z(Ti) = l(T 2 ) and that Ti and T 2 don't coincide in a neighborhood of g* due to the fact that 
they are two different optimal geodesies. Now consider the two projections T\ and £2: they are trajectories 
in L(4,l) steering [Id] to [g] and their length is l(Tj) = d(g*,Xi) = Tmn x e ^d(g*,Xj) due to the fact that 
d(g*,xi) = d(g*,x 2 ) < d(g*,x*). Hence Tj_ and £2 are optimal geodesies. Notice that they are distinct in a 
neighborhood of [g] because they are projections of two trajectories that are distinct in a neighborhood of g*. 
Observe that the choice of g* G [g] is not relevant, as a consequence of Theorem [6] 

We end proving that there is no cut point outside K loc U K sym . By contradiction, let [g] G L(4, l)\(K loc U 
K syrn U {[Id]}) be a cut point: thus there are 71 and 72 two different optimal geodesies of L(4, 1) steering [Id] 
to [g\. Consider their two lifts 71 and 72 starting from Id G SU(2). They reach respectively gi,g 2 £ [g] and 
they are optimal in 577(2). We consider two distinct cases: 

• 9i — 92'- we have two distinct optimal trajectories 71 and 72 steering Id to the point g\. Hence g\ G if id- 
Observe that g\ [Id] otherwise [g] = [Id]. Thus g\ G ifid\[Id], hence [g] G if^ oc . Contradiction. 

• 9i 7^ 92 • consider the two lifts of 71 and 72 reaching both the point g\ and call them ji and 72 respectively. 
Observe that 71 = 71, hence 71 (0) = Id. Instead, 72(0) = x 2 7^ Id: we prove it by contradiction. If 
72(0) = Id, then 72 = 72 by the uniqueness of the lift. Hence 72 (T) = g 2 7^ gi, that contradicts the 
definition of 72. Contradiction. Observe that x 2 G [Id] due to the fact that 72(0) = [Id]. 

We have that both 71 and 72 are optimal, hence d(Id, #i) = /(71) = K71) = '(72) = K12) — d(x 2l gi) with 
x 2 G [Id] and x 2 7^ Id. Observe that, choosing any x* G [Id], we have d(Id,#i) < d(x*,gi) due to the 
optimality of 71. Thus [g] G K sym . Contradiction. 

□ 

We now show some pictures of the projections on S 2 of optimal geodesies on L(4, 1), i.e. the optimal 
trajectories for our problem. 

In Figure [71 some optimal trajectories with cusps are shown: they meet in the same point with the same 
direction, where they lose optimality. This point with direction lies in the local cut locus if loc . 
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In Figure [8] we show smooth optimal trajectories, that are indeed half-circles: the points where they lose 
optimality lie in the intersection K loc D K a . 

In Figure [9] we show optimal trajectories (both smooth and with cusps) meeting in one of the two points of 
intersections of the three strata K loc D K a D K b . 




Figure 7: Projections of optimal trajectories with cusps. 




Figure 8: Projections of smooth optimal trajectories. 



A Basic results in sub-Riemannian geometry 

In this Appendix we recall some definitions and results about sub-Riemannian geometry, in particular about 
k p manifolds. For a deeper presentation of sub-Riemannian geometry see e.g. [5} fT5], 
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A.l Sub-Riemannian manifold 

A (n, m)-sub-Riemannian manifold is a triple (M, A,g), where 

• M is a connected smooth manifold of dimension n; 

• A is a Lie bracket generating smooth distribution of constant rank m < n, i.e. A is a smooth map that 
associates to q £ M a m-dim subspace A(g) of X^M, and V £ M we have 

Bpan{[/i, [. . . [fk-ufk] ■ ■ •]](<?) I /; € Vec(M) and /;(p) e A(p) VpeM} = T g M. (9) 
Here Vec(M) denotes the set of smooth vector fields on M . 

• g q is a Riemannian metric on A(g), that is smooth as function of q. 

The Lie bracket generating condition (|9]) is also known as Hormander condition. 

A Lipschitz continuous curve 7 : [0, T] — > M is said to be horizontal if 7(f) £ A(j(t)) for almost every 
t £ [0,T]. Given an horizontal curve 7 : [0, T] — > M, the length 0/7 is 

Z W = / vg^rni^rft. (io) 

The distance induced by the sub-Riemannian structure on M is the function 

^(<7o, 9i ) = inf {K7) I 7(0) = 7( T ) = Qi, 7 horizontal}. (11) 

The hypothesis of connectedness of M and the Lie bracket generating assumption for the distribution 
guarantee the finiteness and the continuity of •) with respect to the topology of M (Chow's Theorem, see 
for instance [4]). The function ■) is called the Carnot-Charateodory distance and gives to M the structure 
of metric space (see [5] fl5]). 

It is a standard fact that £(7) is invariant under reparameterization of the curve 7. Moreover, if an admissible 
curve 7 minimizes the so-called energy functional 

£(7)= / &y(t)(7W,7(*))*. 
Jo 
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with T fixed (and fixed initial and final point), then v = \/S7(t)(7(£)> 7W) * s constant and 7 is also a minimizer 
of /(•). On the other side a minimizer 7 of /(•) such that v is constant is a minimizer of E(-) with T = l{p/)/v. 

A geodesic for the sub-Riemannian manifold is a curve 7 : [0, T] — > M such that for every sufficiently small 
interval [£1,^2] C [0, T], 7| [t t j is a minimizer of £"(•)• A geodesic for which g 7 (t) (7(^)5 7(0) is (constantly) 
equal to one is said to be parameterized by arclength. 

Locally, the pair (A, g) can be given by assigning a set of m smooth vector fields that are orthonormal for 
g, i.e. 

A(q) = span {F 1 (q), . . . , F m (q)} , g q (F z (q), F 3 (q)) = 5 ZJ . (12) 

When (A, g) can be defined as in ([T2]) by m vector fields defined globally, we say that the sub-Riemannian 
manifold is trivializable. 

Given a (n, m)- trivializable sub-Riemannian manifold, the problem of finding a curve minimizing the energy 
between two fixed points qo,q\ G M is naturally formulated as the optimal control problem 

m „t m 

q = y2uiFi(q), Ui ER, / V^ 2 (£)^^min, q(0) = q , q(T) = Ql . (13) 
i=i Jo i=i 

It is a standard fact that this optimal control problem is equivalent to the minimum time problem with 
controls ui, . . . , Um satisfying u\ + • • • + < 1. 

When the manifold is analytic and the orthonormal frame can be assigned through m analytic vector fields, 
we say that the sub-Riemannian manifold is analytic. 

The k p manifold presented below are trivializable and analytic since they are given in terms of left- 
invariant vector fields on Lie groups. 

A. 2 First order necessary conditions, Cut locus, Conjugate locus 

Consider a trivializable (n, m)-sub-Riemannian manifold. Solutions to the optimal control problem (fT3|) are 
computed via the Pontryagin Maximum Principle (PMP for short, see for instance [H [H [161 [21] ) that is a first 
order necessary condition for optimality and generalizes the WeierstraB conditions of Calculus of Variations. 
For each optimal curve, the PMP provides a lift to the cotangent bundle that is a solution to a suitable pseudo- 
Hamiltonian system. 



Theorem 5 (Pontryagin Maximum Principle for the problem ([T3])). Let M be a n- dimensional smooth manifold 
and consider the minimization problem / TlffJ) . in the class of Lips chit z continuous curves, where Fi, i = 1, . . . , m 
are smooth vector fields on M and the final time T is fixed. Consider the map : T*M x R x R m — > R defined 
by 

m m 

Jif(q,X,po,u) := < \ ^UiFi(q) > +p ^2u 2 i (t). 

1=1 2=1 

If the curve q(.) : [0,T] — > M corresponding to the control u(.) : [0,T] — > R m is optimal then there exist a never 
vanishing Lipschitz continuous covector A(.) :t e [0,T] h \(t) g T*^M and a constant po < such that, for 
a.e. t G [0,T]; 

(i) = -^-((zW,A(t),po,tz(t)), 

(ii) \( t )=-^(q(t) : \(t) :Po Mt)), 

(iii) ^(g(t),A(t),p ,^)) = 0. 
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Remark 6. A curve q(.) : [0,T] — > M satisfying the PMP is said to be an extremal In general, an extremal 
may correspond to more than one pair (A(.),po)- If an extremal satisfies the PMP with po 7^ 0, then it is called 
a normal extremals. If it satisfies the PMP with po = it is called an abnormal extremal. An extremal can be 
both normal and abnormal. For normal extremals one can normalize po = —1/2. 

If an extremal satisfies the PMP only with po = 0, then it is called a strict abnormal extremal. If a strict 
abnormal extremal is optimal, then it is called a strict abnormal minimizer. For a deep analysis of abnormal 
extremals in sub-Riemannian geometry, see [6j [12] . 

It is well known that all normal extremals are geodesies (see for instance [4]). Moreover if there are no strict 
abnormal minimizers then all geodesies are normal extremals for some fixed final time T. This is the case for 
the so called 3-D contact case, i.e. a sub-Riemannian manifold of dimension 3 for which A = span {F±, F2} 
where F±, F2 is a pair of vector fields such that for all q G M, span {Fi(q), i 7 ^), [Fi(q), i^Oz)]} = T q M. Indeed 
for contact structures there are no abnormal extremals (even non strict). 

In this case, from (iii) one gets Ui{t) =< A(£),i^(£) >, i = 1 . . . ,ra and the PMP becomes much simpler: 
a curve q(.) is a geodesic if and only if it is the projection on M of a solution (\(t),q(t)) for the Hamiltonian 
system on T*M corresponding to 

1 m 

H(\q) = -(J2 < \Fi(q) > 2 ), q € M, A G T*M. 

1=1 

satisfying £T(A(0), q(0)) ^ 0. 

Remark 7. Notice that H is constant along any given solution of the Hamiltonian system. Moreover, H = ^ if 
and only if the geodesic is parameterized by arclength. In the following, for simplicity of notation, we assume 
that all geodesies are defined for t G [0, +00). 

Fix q G M. For every A G T* Q M satisfying 

ff(A ,<Zo) = l/2 (14) 

and every t > define the exponential map Exp(Ao, t) as the projection on M of the solution, evaluated at time 
t, of the Hamiltonian system associated with iJ, with initial condition A(0) = Ao and g(0) = qo. Notice that 
condition (HD defines a hypercylinder A qo - S 771 ' 1 x R n - m in T* Q M. 

Definition 2. The conjugate locus from qo is the set C qo of critical values of the map 

Exp- A ^ XR+ M > 

P ' (A ,t) ^ Exp(A ,t). 

For every Ao G A qo , let t(\o) be the n-th positive time, if it exists, for which the map (Ao,t) Exp(Ao,£) is 
singular at (Ao,t(Ao)). The n-th conjugate locus from qo C^ Q is the set {Exp(Ao, t(Ao)) | t(Ao) exists}. 
The cut locus from qo is the set K qo of points reached optimally by more than one geodesic, i.e., the set 

K, . {, e M 1 1 A„ A, e K , A, ^A 2 , t e «* «* ^'^W » HU } 

Remark 8. It is a standard fact that for every Ao satisfying (fT4|) . the set T(Ao) = {? > | Exp(A, t) is singular 
at (Ao,f)} is a discrete set (see for instance [3]). 

Remark 9. Let (M, A, g) be a sub-Riemannian manifold. Fix qo G M and assume: (i) each point of M is 
reached by an optimal geodesic starting from qo] (ii) there are no abnormal minimizers. The following facts 
are well known (a proof in the 3-D contact case can be found in [3]). 

• the first conjugate locus Cq Q is the set of points where the geodesies starting from go lose local optimality; 

• if g(.) is a geodesic starting from qo and i is the first positive time such that q(t) G K qo U C^ o , then q(.) 
loses optimality in f, i.e. it is optimal in [0, t\ and not optimal in [0, t] for any t > i; 
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• if a geodesic q(.) starting from qo loses optimality at t > 0, then q(t) G if go U G* o ; 

As a consequence, when the first conjugate locus is included in the cut locus, the cut locus is the set of points 
where the geodesies lose optimality. 

Remark 10. It is well known that, while in Riemannian geometry K qo is never adjacent to go, in sub-Riemannian 
geometry this is always the case. See [2]. 

A. 3 k © p sub-Riemannian manifolds 

Let L be a simple Lie algebra and Kil(X, Y) = Tr(adx ° ady) its Killing form. Recall that the Killing form 
defines a non-degenerate pseudo scalar product on L. In the following we recall what we mean by a Cartan 
decomposition of L. 

Definition 3. A Cartan decomposition of a simple Lie algebra L is any decomposition of the form: 



Definition 4. Let G be a simple Lie group with Lie algebra L. Let L = k® p be a Cartan decomposition of~L. 
In the case in which G is noncompact assume that k is the maximal compact subalgebra of~L. 

On G, consider the distribution A(g) = gp endowed with the Riemannian metric g^(^i, V2) =< 9~ lv i-> 9~ lv 2 > 
where < , >:= a Kil|p( , ) and a < (resp. a > 0) if G is compact (resp. non compact). 

In this case we say that (G, A, g) is a k ® p sub-Riemannian manifold. 

The constant a is clearly not relevant. It is chosen just to obtain good normalizations. 

Remark 11. In the compact (resp. noncompact) case the fact that g is positive definite on A is guaranteed by 
the requirement a < (resp. by the requirements a > and k maximal compact subalgebra). 

Let {Xj} be an orthonormal frame for the subspace p C L, with respect to the metric defined in Definition 
IU Then the problem of finding the minimal energy between the identity and a point g\ G G in fixed time T 
becomes the left- invariant optimal control problem 



Remark 12. This problem admits a solution, see for instance Chapter 5 of [TP] . 

For k © p sub-Riemannian manifolds, one can prove that strict abnormal extremals are never optimal, 
since the Goh condition (see [3]) is never satisfied. Moreover, the Hamiltonian system given by the Pontryagin 
Maximum Principle is integrable and the explicit expression of geodesies starting from go and parameterized 
by arclength is 



where Ak G k, A p G p and < A P ,A P >= 1. This formula is known from long time in the community. It was 
used independently by Agrachev [1], Brockett [TT] and Kupka (oral communication). The first complete proof 
was written by Jurdjevic in [17] . The proof that strict abnormal extremals are never optimal was first written 
in 0. See also [HCEE]. 

Acknowledgements: We thank A. Agrachev for his help in recognizing the structure of the projective 
tangent bundle. We thank L. Paoluzzi for many explanations on lens spaces. We thank M. Saponi and G. 
Senaldi for their help with pictures. 



L = k © p, where [k, k] C k, [p, p_ 



i]Ck, [k,p]Cp. 




Uj G L°°(0,T) , / 5>^) dt - min, g(0) = Id, g(T) = 9l . 




g(t) = g e^ A ^e- 
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